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Abstract 

Spacecraft attitude control using only magnetic torques is a time-varying system. Many designs 
were proposed using LQR and Hoo formulations. The existence of the solutions depends on the 
controllability of the linear time-varying systems which has not been established. In this paper, we 
will derive the conditions of the controllability for this linear time-varying systems. 

Keywords: Spacecraft attitude control, linear time-varying system, reduced quaternion model, con¬ 
trollability. 


Office of Research, NRC, 21 Church Street, Rockville, 20850. Email: yaguang.yang@verizon.net 



1 Introduction 


Spacecraft attitude control using magnetic torque is a very attractive technique because the implementa¬ 
tion is simple, the system is reliable (without moving mechanical parts), the torque coils are inexpensive, 
and their weights are light. The main issue of using only magnetic torques to control the attitude is 
that the magnetic torques generated by magnetic coils are not available in all desired axes at any time 
PP. However, because of the constant change of the Earth’s magnetic held as a spacecraft circles around 
the earth, the controllable subspace changes all the time, many researchers believe that spacecraft’s at¬ 
titude is actually controllable by using only magnetic torques. Numerous spacecraft attitude control 
designs were proposed in the last twenty five years exploring the features of the time-varying systems 
[a Ei mEusuzusuni unum na m] . Some of these papers tried Euler angle model and Linear Quadratic 
Regulator (LQR) formulations [2j 31 [4, 0 111 ] which are explicitly or implicitly assumed that the con¬ 
trollability for the linear time-varying system holds so that the optimal solutions exist M- But for the 
problem of spacecraft attitude control using only magnetic torque, no controllability condition has been 
established for this linear time-varying system to the best of our knowledge. 

Other researchers [SE 1 proposed direct design methods using Lyapunov stabilization theory. The 
existence of the solutions for these methods implicitly depends on the controllability for the nonlinear 
time-varying system. Therefore, Bhat m investigated controllability of the nonlinear time-varying sys¬ 
tems. However, the condition for the controllability of the nonlinear time-varying systems established in 
this paper is hard to be verified and is a sufficient condition. 

Recently, a reduced quaternion model was proposed in [16] and its merits over Euler angle model were 
discussed in mmm- The reduced quaternion model was also used for the design of spacecraft attitude 
control system using magnetic torque nnnuna. Because the controllability of the linear time-varying 
systems was not established, the existence of the solutions was not guaranteed. 

In this paper, we will consider the reduced linear quaternion model proposed in |16| and incorporate 
control model using only magnetic torques. We will establish the conditions of the controllability for this 
very general linear time-varying system. The same strategy can easily be used to prove the controllability 
of the Euler angle based linear time-varying system considered in [5]. However, we will not derive the 
similar result because of the merits of the reduced quaternion model as discussed in [H HO [IB]. 

The remainder of the paper is organized as follows. Section 2 provides a description of the linear 
time-varying model of the spacecraft attitude control system using only magnetic torque. Section 3 gives 
the proof of the controllability for this linear time-varying system. The conclusions are summarized in 
Section 4. 


2 The linear time-varying model 

Let J be the inertia matrix of a spacecraft defined by 
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We will consider the nadir pointing spacecraft. Therefore, the attitude of the spacecraft is represented 
by the rotation of the spacecraft body frame relative to the local vertical and local horizontal (LVLH) 
frame. Therefore, we will represent the quaternion and spacecraft body rate in terms of the rotations of the 
spacecraft body frame relative to the LVLH frame. Let u = [uj\, u> 2 , W 3 P be the body rate with respect to 
the LVLH frame represented in the body frame, loq be the orbit (and LVLH frame) rate with respect to the 
inertial frame, represented in the LVLH frame. Let q = [go, qi, <Z 2 , <Z 3 ] T = [go,q T ] T = [cos(§), e T sin(^)] T 
be the quaternion representing the rotation of the body frame relative to the LVLH frame, where e is 
the unit length rotational axis and a is the rotation angle about e. Therefore, the reduced kinematics 
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Assume that the inertia matrix of the spacecraft is diagonal which is approximately correct for real 
systems, let the control torque vector be u = [u x ,u y ,u z ]' T , then the linearized nadir pointing spacecraft 
model with gravity gradient disturbance torque is given as follows m : 
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where 


fa = [8(J33 — J22)uq\/J u ( 4 a) 

/46 = (—Ju + J 22 — ^ 33 )^ 0 /J 11 (4b) 

/64 = {Jll — J22 + ^33) w o/^33 ( 4 c) 

fb2 = [6( J 33 - Jii)uq]/ J22 (4d) 

/63 = [ 2 (Jn — ^22)^0]/J33. ( 4 e) 


The control torques generated by magnetic coils interacting with the Earth’s magnetic held is given by 
(see [1]) 

u = m x b 

where the Earth’s magnetic held in spacecraft coordinates, b(f) = [61(f), 62(f), 63(f)] T , is computed using 
the spacecraft position, the spacecraft attitude, and a spherical harmonic model of the Earth’s magnetic 
held [l 9 j : and m = [mi, m2, m3] 1, is the spacecraft magnetic coils’ induced magnetic moment in the 
spacecraft coordinates. 

The time-variation of the system is an approximate periodic function of b(t) = b(t + T) where T = 
is the orbital period. This magnetic held b(t) can be approximately expressed as follows j 5 j: 
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where i m is the inclination of the spacecraft orbit with respect to the magnetic equator, [if = 7.9 x 10 15 
Wb-m is the held’s dipole strength, and a is the orbit’s semi-major axis. The time f = 0 is measured at the 
ascending-node crossing of the magnetic equator. Therefore, the reduced quaternion linear time-varying 
system is as follows: 
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Substituting © into © yields 


where 


Therefore, we have 
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In matrix format, we have 
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A special case is when i m = 0, i.e., the spacecraft orbit is on the equator plane of the Earth’s magnetic 
field. In this case, h(t) = [0, — ^r,0] T is a constant vector. The linear time-varying system of this special 
case is reduced to a linear time-invariant system whose model is given by 


9i 


0 

0 

0 

.5 

92 


0 

0 

0 

0 

93 


0 

0 

0 

0 

U 1 


fil 

0 

0 

0 

U>2 


0 

fs2 

0 

0 

. 


0 

0 

f&3 

/64 


= Ax + Bm. 


0 

0 


9i 


0 

.5 

0 


92 


0 

0 

.5 


93 

+ 

0 

0 

/*46 


UJi 

0 

0 

0 


0J2 


0 

0 

0 


. W 3 


b 2 / J 33 


0 

0 

0 

0 

0 

0 


0 

0 

0 

—b 2 /Jn 
0 
0 


m i 
to 2 
to 3 


(13) 


3 The proof of the controllability 


The definition of controllability of linear time-varying systems can be found in m page 124]. 

Definition 3.1 The linear state equation m is called controllable on [to,tf] if given any Xq, there exists 
a continuous input signal m (t) defined on [to,tf] such that the correspnding solution of ( 0 ) satisfies 
x{t f ) = 0 . 

A main theorem used to prove the controllability of (j 6 j) is also given in [2Ql page 127]. 

Theorem 3.1 Let the state transition matrix 4>(t, t) = e A ( t_T b Denote 


K,Ct) = [<!>(*, t)B(t)] 


3 = 1,2,. 


(14) 


if p is a positive integer such that, for t € [to,tf], B(i) is p time continuously differentiable. Then, the 
linear time-varying equation © is controllable on [t-o,tf] if for some t c € [to,tf] 


rank[K 0 (t c ),Ki(t c ),... , K p (f c )] = n. 


(15) 


Remark 3.1 If A and B are constant matrices, the rank condition of ESP for the linear time-varying 
system is reduced to the rank condition for the linear time-invariant system !J?0i page 128], i.e., if 

rank [B, AB,.... A" -1 B] = n. (16) 

then the linear time-invariant system (A, B) is controllable. 


First, we consider the special case of m the time-invariant system when the spacecraft orbit is on the 
equator plane of the Earth’s magnetic field ( i m = 0). Let X denote any 3x3 diagonal or anti-diagonal 
matrix with the second row composed of zeros 
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and A denote any 3x3 diagonal matrix with the form 
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It is easy to verify that if Xj £ X, Xj £ X, and A/. £ A, then X,Xj £ X, X; + Xj £ X, and A^X; £ X. 
Using this fact to expand the matrix [B, A 2 B, A 3 B, A 4 B, A 5 B], where A and B are defined in m, 
shows that the second row of the controllability matrix in (fl 6 l) is composed of all zeros. This proves that 



















if the spacecraft orbit is on the equator plane of the Earth’s magnetic field, the spacecraft attitude cannot 
be stabilized by using only magnetic torques. 

Now we show that under some simple conditions, the linear time-varying system © is controllable 
for any orbit which is not on the equator plane of the Earth’s magnetic field, i.e., i m ^ 0. From (TUEl) . we 
have 

K 0 (t) = S(t,t)B(t) = e A ^- t} B(t) = B(t), 
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To show that this matrix is full rank for some t c , we show that there is a 6 x 6 submatrix whose determinant 
is not zero for wo t c = f • In view of ©, ©, and ifTUI) . for this 6 C , we have 
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Therefore, in view of Theorem 13.11 the time-varying system is controllable if 


/64^>42(ic) - 26g 2 7^ 0, 


and 


^51^62/46^61 — ^42 ^-/52&51 + 651^ ^61 + 2/63^61^42^51 7^ 0. 
Using ©, ©, ©, m, and noticing that sin(u;cJc) = sin( : |) = 1 , we have 


J33 

2/4/icq sin(i m ) 
a 3 ( Jn J 33 ) 

the first condition C5l) is reduced to 


/64^42(ic) — 26g 2 

W, - -fe + JssVo 2 W sin(im) _ 2 i^i sill(jm) 


a 3 Jn 

(•/s3 ~ ^ 22 ), 


a 3 J; 


33 


J 33 7 ^ <^22 • 


Repeatedly using the same relations, we have 


/ / 2/if 

^51^62/46^61 = ( ~ q 3j sm ' 

,3, ,2 


/• 0 ■ f (-Jll + J22 — J33 )uq\ 

M ) lSfa sm( ’" ) J l-fc-j 


2/ijO/ 0 (—Jn + J 22 — ^ 33 ) 
a 9 J 11 J 22 J33 


sin (im)eos(zm), 


( 25 ) 

( 26 ) 


( 27 ) 


- COs(im) 
a 3 J 33 


( 28 ) 


— &42 ( ^/SS^SI + 5 51 ) ^61 — 


2 /i/ . A / 3 (J33 


— Jn)wg 


2/X/ sin(* m ) N ) + 2/ ^ /a; ° sin(z„ 


a 3 J; 


a 3 J 2 


M/ 

a 3 J; 


33 


cos(i m ) 


2/i/ sin(z m )^ { ^TfT sin(i TO )(-3J 33 + 3Jn + J 22 ) 


a 3 Jn 

a 3 ^33 

,3, ,2 


cos(z m ) 


4/i/W 0 (— 3 J 33 + 3Jn + J 22 ) . 2 


a 9 Jn J| 2 J 33 


sin (i m ) cos(i m ), 


and 


1 


f63^61^32^31 = 


J ‘ 22 ^ >bjQ f -cos(z m )^) sin(i m ) N j f — (i m ) 

J33 V a 3 J 33 / Va 3 Jn / V a 3 J 22 


4ft 3 WQ(Jn — J 22 ) . 


a 9 Jn J 22 J 33 

Combining (Ell), (EU), and (®|), we can rewrite 


sin (z m ) cos(i m ). 


as 


^51^62/46561 — ^42 ^-/52551 + &61 + 63 ^ 61 ^ 32^51 

u 3 w 2 

g T ?2 t 2 “ S t n (*m) cos (i m ) 
a ^11^22^33 

(2 J 22 ( —Jll + J 22 — J 33 ) + 4J33( —3 J 33 + 3 Jll + J 22 ) + 4J 2 2( Jll — J 22 )) 

2 /i?WQ „ 

— 3 —- 2 2 sin (i m ) COs(i m )[J 22 (J 11 — J 22 + ^ 33 ) ~ 6^3^33 — >/ll)]- 

a <711^22^33 


( 29 ) 


(30) 


(31) 


























Therefore, the second condition of (EE1) is reduced to 


•^ 22(^11 — J22 + J33) 7 = 6 J 33 (•/as — Jn)- (32) 

We summarize our main result of this paper as the main theorem of this paper. 

Theorem 3.2 For the linear time-varying spacecraft attitude control system 0) using only magnetic 
torques, if the orbit is on the equator plane of the Earth’s magnetic field, then the spacecraft attitude 
is not fully controllable. If the orbit is not on the equator plane of the Earth’s magnetic field, and the 
following two conditions hold: 


J33 i 1 J22, (33a) 

^ 22(^11 — J22 + J33) 7 ^ 6 ^ 33(^33 ~ Ji 1 ), (33b) 

then the spacecraft attitude is fullly controllable by magnetic coils. 

Remark 3.2 The controllability conditions include only the spacecraft orbit plane and the spacecraft 
inertia matrix which are very easy to be verified. 

4 Conclusions 

In this paper, the controllability of spacecraft attitude control system using only magnetic torques is 
considered. The conditions of the controllability are derived. These conditions are easier to be verifed 
than previously established conditions. 
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